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Ali~it rail . A tlilrd order, nonlinear, singularly perturbed oiitimal control 
probletn Is considered under assumptions whlcli assure that the full problem 
Is slniiular and tlie reduced problem is nonsingular. The separation between 
the singular arc of the full problem and the optimal control law of the 
leduced one, both of which are hypersurfaces In state space, Is of the same 
order as the small parameter of the problem, Boundary layer solutions are 
constructed which are stable and roach the outer solution in a finite Ime, 
A uniformly valid composite solution is then formed from the reduced and 
boundary layer solutions. The value of the approximate solution Is that it 
Is relatively ease to obtain and does not Involve singular arcs. To Illus- 
trate the utllltv of ttie results, the technique Is used to obtain an 
approximate solution of a simplified version of the aircraft minimum time- 
to-cllmb probh A numerical example is Included. 

Keywords . Singular perturbations, optimal control, singular arcs, flight 
mechanics. 


INTRODUCTION 


Singular arcs h.ivc been found to occur In 
optimal control problems from many fields. 

They are particularly prevalent In flight 
meclianlcs because in many problem formula- 
tions vehicle thrust Is mode-led as a bounded 
control variable which appears linearly in 
the state equations and performance Index. 
Well-knawn examples are the problems of maxi- 
mizing the range of a rocket in vacuo 
(Leltmann, I96b), maximizing the range of a 
lifting rocket in horizontal flight (Leltmann, 
1966), minimizing the fuel required for orbit 
rendezvous of a rocket (Breakwell and Dixon, 
1975), Goddard's problem of maximizing the 
final altitude of a sounding rocket (Bryson 
and Ho, 1969), and minimizing the tlme-to- 
cllrab of an airplane under certain approxi- 
mations (Ardema, 1975). Current Interest in 
flight mechanics problems with possible 
singular arcs is evidenced by many recent 
papers concerning the optimality of cruise 
for airplaies (Schultz, 1974; Speyer, 1976). 

In this paper, we investigate the usefulness 
of singular perturbation techniques for 
dealing with singular arc problems by 
analyzing a relatively low-order but other- 
wise general system. This system encompasses 
many flight mechanics problems Including 
Goddard's problem and a version of the mini- 
mum tirae-to-climb problem. 


dx 

ill 

dt 


f(x,y,z) i'(x,y,z)u 
g(x,y,a) + g' (x,y,z)u 


c 


dz 

dt 


h(x,y ,z) + h' (x,y ,z)u 


( 1 ) 


subject to sultiible boundary conditions, 
where x, y, and z are scalars, c > 0 Is a 
parameter, and tf Is unspecified. It is 
desired to 

Min ^ l^x(tf) ,y (tf )j (2) 

The scalar control is restricted to 

u^<u<uj5 (3) 

Since the Hamiltonian for this problem will 
be linear In u, a singular arc is possible. 
It is assumed that a unique optimal control 
exists and that In the region of state space 
of interest:' 

A. fxi fy< ^ X' ^ y ’ ^ Z ’ 8x • ky 1 8z • 

R'x’ s'y» s'z» ^x* l*'x> h'y* h'z 

are defined and continuous. 

B. h+h'u„,-!^0, h + h'u(,j»i'0, 

sgn (h + h'u^) = - sgn (h + h'u^). 


Consider the third order system, defined on 'Functional dependence will be omitted when 

0 < t < t£, this docs not result in confusion. 
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I . i.jtiH-r - f' ) * () 

(k - «' II or both. 

o. llio 'iliu’ii ;.ir .iri. ati.stU*?. ( i) .inJ the 
gfiicriil I /(.'d roim-xlix oon>litlon Uifvwon 
and Ho, I9h‘i; Kojip, .jnd Moyer, 

Ihus wx' havv a xlnguiarly porturbod optimal 
.'ontrol pr"|i|> m tor will, ii the J‘iill problem 
1-s HinKiil.ir bat the fodut ed problem 1« not. 

A.s>iiraptlon B .ifihurt'st th.il 1) h' * Q, i) i. ••an 
hi’ eltlier dfi rt’.isfed or Increatied , and 3) the 
alternative procedure may be used for the 
reduced problem tArdema, 197 j). Assumption C 
tliem implies that the reduced problem Is not 
Itself SltiKulur. 

Singular perturbation problems involving sln- 
gulai arcs have been considered previously, for 
example by O'Malley (1975) -nd O'Malley and 
Jameson (1975). However, tliese references 
consider linear systems for which tlie full 
problem is aonsingular but the reduced prob- 
lem is singular; therefore thev bear little 
connection with the present work. Rather, 
our approach follows Calls© (197 it who first 
proposed using singular perturbation tech- 
niques to solve nonlinear singularly per- 
turbed optimal control problems. This was 
given furtlier consideration by Ardema (1975), 


ANALYSIS 

The Hamiltonian for the problem jt.st stated 
Is 

H » (f + f'u)\ + (g + g’u)Ay + (h + li'u)ijj 


and the 

adjoint 

oquat Ions 

are 


dX X 





~dt “ 

-(f^4 

fj/u)Xx - 

^Hx 

gx'u)'y 


- (hx + 

hx'.O’z 



dt 

-dv + 

fy'u)Xx - 

(gy 

gy'u)Xy 


- (\ b 

h'yU)X2 



d'z _ 
" df 

-(f2 + 

f2*u)Xx - 

(gz + 

g',U))y 


(h2 4 

h'2U)'2 



The "switching 

funct ion" 

(coefficient - 

In (4) 

is 





S = 

fXx + g’x 

V + h' 



The part of H not Involving S is 

H = fX^ + gXy + hXj, (7) 

Using (1) and (5), the time derivative of (6) 
is, after multiplying by c, 

= [f,*h - h'f, + e(f^'f + fy’g - ff^ 

- g'fyll^x + (gi:’h - h'g^ -t t(g^'f 

+ gy'g - f'gx “ g'gy)J^y + 

+ c (h'^f + h'yg - f hx - g’hy)1>z (8) 


liie terms Involving u t in. •! . ^ii . on i 
-Singular arc ..i an iiptlmal t r iln tor-. , H - 0 
,ind S “ 0 , HI* th.it We hav< till |.»i lowing 

SVStem of ht>it!’’gi tleOUS eipiat huiti In till- 
adl-'lnt varl ihli s! 

H • (I ^ 


S - I) 



( 9 ) 


lor omp.it ihl 1 1 tv M livit Is, existence ot non- 
trivial solulJoilsf tile determlti-inl "1 the 
coefflilent m.tti ix of (9i must be /.oro whl* h 
gives 

(fg’ - f'g)0''?h - h'll2 + cUi'xf + hy'g 

- f’h^ - g'h^.t] 4 (gh' - hg')[f^'h - h'f^ 

4 C(f„'f 4 fy'g - f'fx - g'fv)] (hf 

- fli')[gj,'h - h’gg + ’ (gg'f Hy'g 

- f'Kx - g'Ky)] • 0 
which Is the singular arc. 


We next sliow directly that the control law 
for the reduced problem is (10) with c set 
to zero. V.'e adopt the alternative procedure 
and set t • 0 in (11) and eliminate u to 
get the reduced system of equations 

, , 1‘0 

dt * o - 0 no' , 

‘IYo , ^0 

dt " " 8o 


where, for example, f^ » 
performance Index becomes 






The 


's' problem has control variable 


7. and Is nofislngular. For optimal control 


»o “ 




K'no 

+ (go - Ho’ *■ “ 



For compatibility of these equations,^ 


^Since for some combinations of performance 
Index and terminal boundary conditions It 
Is possible that Xq 0 even If = Xy 

^ >2 ~ 0 (Leltmann, 1966), the trivial solu- 
tion of (9) may in fact satisfy the Maximum 
Principle. The subsequent development Is 
tacitly restricted to problems for which 
this does not occur, 

comment similar to the one In footnote^ 
also applied here. 
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* * iiHn “ 'o -'o' ■*■ 

• ''o’».^,)‘«u''o' * 

- - i„h.,') - i> d o 

wlilili Ih tlLU wltli >- * 0. Th'ifs tin- uiriKul.jr 
ar.- of tlu full problem nnJ tho lantrol lav 
of i lu- re^iui fj pri'ii I fin ill f (or by .it m>>Nt 
1 1 ruin i‘t .■rder 


llio t raiiHform,it lou <ip)iroaoli (u ?itii>'.ular ar> .s 
Jove loped (n Kolhv, kupp, and Mu.'or (I9h7) 

Is part lou'arly .ulvantagoouB f«>r slngul.irlv 
porturbed problems. Coder f.;irlv nil Id 
asBUfflpCions, a tran.sfurmat Ion of state vari- 
ables mav be found such that In the new 
varlnblis the state equations are of the form 


dx 

dt 


f (x,y,.'l 


dv 

dt 


« 


glx.v.zi 


dz 

dt 


o(x,y,z) + h' (x,y,z)u 



This transforiuat ion fKelley, Kopp, and Moyer, 
19fc7) may be obtained os follow.si Let the 
solution of 


f * (x.V.z) , ^ g' (x,^,z) 

d2 h'(55,^,2) ’ dz h' (x,y,£) 

be denoted by 
Cx “ 

Cy • ifry(x,y,2) 

where Cx and Cy are constants of integration 
and and i|iy are two independent integrals. 
Then the transformation of variables 

X * 

9 “ !jly(X,y,Z) 
z -• z 

transforms (1) into C15). The advantage of 
this approach for singularly perturbed 
problem.^ is that from (10) and (lA) we See 
chat for (15) the singular arc of the full 
problem and the control law of the reduced 
problem are Identical and are given by 

fgz - gfz * 0 (19) 

The principal difficulty with Implementing 
thlis approach lies in solving the generally 
nonlinear equations (16). 

In many applications, a small parameter must 
be Inserted artificially to create a singu- 
larly perturbed problem (Aiken and Lapidus, 
1976; Ardema, 1975. 1976; Calisv, 1973; and 
Kelley, 1973). In this case, tlie transfor- 
mation approach has the further advantage of 
clarifying where to Insei f the parameter. 

The transformation results in a system in 
which variables x and y are controlled by a 
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n latlvt lv f.ihi I r v.ii table t wliU h Is, In 
turn, • ■•■•t 1 >•! led (ly •! »tllt fu‘.t<r variable 
It, llius it 1 ■ led tbit I If. to bt 

Insi-rtiil •t'« < multiplier of l(u- derlv.Hlvi- 
1. rii. Hi till l.ivt -.t.ite variable . 

llie .fliiive result', are also V.illd for .1 
sllghtlv diiKi.iit problem. Instead ol (1) 
-iiid ( ') .onsIib*r I iu' systi’m 

*[* ‘ 1 (X./1 + 1 ' (x,z)u 
<1 1 

, ‘ J 4. * J / 

Y = h (x,z ) + h * (x ,3t)u 
ui 

j 

wiu-re tl in t>» 

Min ii[x(Lf)] + /q'’*^IL(X,Z) + L'Cx,z)u)dti 

(2>) 


.subject to (3). Then the H.imlltunlan is 

I! * (f + f'ii)'x + (h + h'u)‘^ + 1 1, + l.'u) '^, 

(22) 

Comparing (6) with (22) we see that the same 
results are valid as before provided t., L', 
and Ay are Identified with g, g' and >y, 
respectively .‘* 

Returning to tiie problem defined by (l)-(3), 
the initial zeroth onk-r boundarv laver 
equation (7iiBLI.) Is formed from (1) by trans- 
forming ^ j ** 7^ setting I » U. The 
result is 

dzi , , 

37^ = ^ +h;'ui, (23) 

where, for example, hj *• h(x(0) ,v(0) ,Zj) . 
Similarly for the terminal ZohLK set 

t 2 • (tp - l )/t: to get 

dzo 

37“ “ “^2 - h^'ua (26) 

where, for example, h^ = h(x(tf) ,y(tf) ,z.) , 


CUNSTRUCTION OF APPROXIMATE 
SOLUTION 

We will now construct an approximate solu- 
tion to the problem (l)-(3) under the stated 
assumptions. To be specific, consider the 
case for whlch^ 

1. x(0), y(0), z(0), X(tf) z(tf) are given, 

2, .1 “ <l> [y (tf) ] , 


“*In this case the trivial solution of (9) 
cannot satisfy the Maximum Principle, 
5Tht*‘.e conditions were choi-en only for the 
pun-'Se of discussion; the method is not 
limited to this case. 
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OHIGINAL PA(JB I& 
OF POUR gUAUTV 


t. llii- optimal t-i>n! ro I h< quvni-< 1 m 
•“M, ••SlSi t’lAK* “M». 

Hv aMMiimpt lun (1 m> m<iy bw Molvvtl tor x; 
the Molutlon 1 h written .in x„ • pIRq.v,,). 
lliitl, liileKrnti.* (11) from ( x(U) , v(0) | , to 
x(t(). Thlb glvoa till xi-rotli order outer 
(re<lu<ed) Molullon, denoted x^, (t), V^, (t). 

.ind ZpU^, i«i will a*i an e»tim.)te oi tf, 
denoted *tjj» •*'>‘1 of v(lf), denoted 
Set'ond, xolve till boundary laver eguatlonx. 
The Initial ZOhl.F. lx 

dx , 

— — • li , ♦ h , ' u„ (25) 

dij I I M 

I'lilx lx liitektrated until the reduced Holut ion 
ia reached: 


A. ‘^^l 

) h(x(0> ,y(0) ,x j ) A h* (x(0) ,y(0) ,X| )u.. 

(26) 


At Tj ■ tj*. u lx xet to -iij/hj’ which 
hoIdN Zj at z^(0). A simliar procedure in 
used for the terminal ZOBLE. This construc- 
tion obviously produces stable boundary layer 
solutions which not only approach the outer 
solution as t ^ ® but reach it in finite 
time. * 

A possible representation of the exact solu- 
tion may be obtained by "patching" the 
unsteady portion of the boundary layer solu- 
tions onto the ends of the outer solution. 

The trajectory time is estimated as 

tf ■ tf + T," c + T,*c by this solution 
‘ p ' o * 2 

and the perfc.'mance Index as ^(y„(tf )]. 

U » Q 

All state variables generally will have 
points of non-dif ferentiablllcy. 

A more satisfying approximate solution is 
available from the method of matched asymp- 
totic expansions. cn this method, "compos- 
ite" solutions are constructed from the outer 
and boundary layer solutions (Ardema, 1975, 
1976; Van Dyke, 1964). The most common of 
these Is the additive composite, the zeroth 
order of which takes the following form; 

X^(t) • Xg(t) 

ya(t) “ Yo<t) /£ 

x^ft) - z^(t) + Zi(J) + ‘ 

- *o(0) - 

Such a solution will satisfy all boundary 
conditions. Because "time stands still" in 
the boundary layer for the slow variables x 
and y, the** are no boundary layer correc- 
tions to thesi. variables to zeroth order. 

The slow variables are approximated by Xj,(t) 
and the tra]ectory time bv tr , the 

final value of y(t) by v(^fQj 
performance index by ^ ['’-('Ol . Thus, If 



I knowlidye ut x(t) li o( ii<' p.iriliulir 
(nten-Ht, till' biiiiiid ir I i\i r i '|u it ui i d 
not he solved. 


Ihe i-ASi t , piit'lied, .iild /iToth older .iddillvt 
■ omposite ■•ulut il'iiM (oi .) iMiii.tl sltu.ition 
.ire shown In Kig. 1. It is cIimi from the 
I onxi rue t ii n lh.it for suf f li lent ly sm.ill 
the latter 1m generalll a un I f urml y-va 1 id , 
/erot h-ord-.-r approxlir-at ion , th.it 1m, that 

xlt)»x^(l)*r(i) I 

yit) - v„(t )■*/(•) I (21!) 

z(t ) • Xj(i ) 1 Ml ) I 

on ■'•'t<t(^^. Higher order .ipproxlroat ions 

may be obt. lined by exp.indlng all variables 
in asymptotic power serijs, matching outer 
and boundary layer ixpansloDH, and forming 
composite solutions (Ardeoia, 1975, 1976; 

Van Dyke, 1964). 

A procedure for solving the exact system 
would be as folli'ws; F.qu.ition., (1) are 
integrated from x(0), y(0), z(0), with 
u ■ u^, until Che slngulai arc (10) Is 
reached. The singular arc portion of the 
trajectory la then obtained by integrating 
(1) subject to (10), At some time t', the 
control is set to um and (1) is integrated 
until x(tf) is reached. If the value of r. 
is z(tf) when tlie value of x is x(tf), 
then the solution has been obtained; if it 
isn't, repetitive solutions of the terminal 
arc must be made by varying t', the time of 
departure from the singular arc, until both 
the terminal conditions are met. Since this 
Iteration is not required for the zeroth 
order additive composite solution, obtaining 
such an approximation requires considerably 
less computational effort than does obtain- 
ing the exact solution. 


MINIMIM TIMF.-TO-CLIMB PROBl.EM 


A well-known problem of flight mechanics is 
that of determining the flight patli which 
gives the minimum time-to-climb between a 
given speed and altitude and another given 
speed and altitude for an aircraft in the 
atmosphere. In one simplified version of 
the problem, the two-state approxisiat ion, 
the rotational dynamics are neglected. This 
leads to the state equations 


^ • v(e,h)F(e,h) 
c ^ - v(e,h)u 


(29) 


where 


u • sin y 

v(e,h) “ /2(l^e - h) 

F(e,h) - ^[T(e,h) - Do(e,h) 

- DL(e,h,L;|L 


> OO) 



/ 
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.Illkl 




t ■ link- 
H • .|l t 1 1 udf* 
n • 

V • v» li»i Uy 
I ■ thriiMt 

I\, • ju-r»» I Hi Jrati 

1*1 • ilr.iK i*’ IHt 

i; ■ Kr-ivlty torn i>t>r unit maaa 

V • ti-iu'rio (••^r unit wi ittlti 
y ■ tllRht pJtli .riiRlr 

1. - lift 

In tills nk'Jcl , m nnil t> art- const. mtw. Fur 
tlu ditr iv<it Ion of tliotf t«c(u.it Iona ,ind .i com- 
plc-tv discussion of thr assumpi ions tli«y 
Imply, the rc.)di-r Is referre-d tu Chapter 2 u( 
Ardi n.i (107''), It Is doslrtd to find the 
control u which gives 

Min jt 

^ 0 

suh)ect to - 1 ^ u < 1 and 

i'(0) • F ,hlO) • H .e(l^) - E<,h(tf) - Hf 

(11) 

The problem Just stated is of the type (20)- 
'21) with 

X • e, r ■ h, f ■ vf, f • 0, h ■ 0, h' ■ V 
I - 0. g - L - 1. g' • L' • 0 

(32) 

All four assumptions listed in the Introduc- 
tion will be satisfied for physically meaning- 
ful problems. The problem is already in the 
form (1)) and c has been inserted in front 
of the derivative of the fast variable h, 

F rom (19), 


. 0 

3h 


(33) 


which is both the singular arc of the full 
ptoblem and the control law of the reduced 
problem. The reduced problem is 

de 

- v(fc„,ho)F(eo,ho) (34) 

subject to 6^(0) • and * ^f where 
Uq “ 0 and h^ is the control, determined from 
(33). The initial and terminal ZOBLE's are 


h ( 


' ■ 'l 



♦ ' ' I. ' 

♦ .F, - 



( <#i) 


In pat : l> ul.ir , the time- to ,>t to the out-r 
■«t>liitl->n .iri- glvt-n h 


I,* [•)„ - hold) - .v„ - ii„] 

I / [.i;, h„(rf) - •T7-‘‘”f] 


(17) 


M.e zeroth order additive compositi approxi- 
mate solution (27) mav be written a- 


e^(t) 

h^(t) 


h„(t) e (-,* 

* j(t ♦ t • - 


- t )Ih, (t )- h (0) 
I o 



>(3H) 


for this problem, where a(t) is the unit 
step function and c has been set to 1, its 
proper value. 


It 1s interesting to compare this approximate 
solution with solutions obtained from differ- 
ent dynamic models. One advantage of the 
singular perturbation viewpoint i» that it 
provides a convenient means for making such a 
comparison. Many of the approximations rom- 
owinly used in flight mechanics art critcally 
discussed in Leltmann (1962). The most 
widely used approximate formulation of the 
minimum t Ime-to-cl Imh problem is that of 
energy-state (Ardema, 1973; Bryson, Desal 
and Hoffman, 1969), In this model, there is 
only one state variable, e; the state equa- 
tion is the first of (29), 


?- • v(e,h)F(e,h) 


(39) 


subject to e(0) “ E„ and e(tf) “ F, where 
h is the control variable. Th « . -si 
control law is 


0 


(40) 


dh) ' 

d77 " v(Eo.hj)u, 

dh, 

— - -v(Ej.h2)u^ 

I 

Assuming that the reduced solution 
"right-hand side" of both boundary 
we have Uj«-1 and - 1 . Using 
are then easily Integrated to give 


^Not to be confused with the function h of 
the preceding equations. 


The path defined by (40) is often called the 
"energy climb path" for obvious reasons. 
Boundary conditions on h are met bv adjoin- 
ing constant energy arcs (traversed in seio 
time in this approximation) to the solution of 
(39) with (40). Because 1) the outer problem 
of the two-state approximation 1s Identical 
to (39) and (40), and 2) boundary layer solu- 
tions (3b) are constant energy arcs, the path 
in (e, h) space obtained by "patching" 
together the outer .and boundary layer solu- 
tions of the two-state approximation is the 
same as that obtained from the energy-state 
approximation. The only difference is that 
the constant energy arcs are not traversed in 
zero time in the patched two-state solution. 


(35) 


lies on the 
condit ions, 
(30), (35) 
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()n tilt* otiii-r tiMml, thi* t lint*a>to-i I leib, vatl- 
nkitfil hv tlif «• roth nrdfr .lUdltlvi cniBpoalti* 
aoiutlon ol tliv two-»t<itt* mtidi'l and bv tbr 
Noliitlun of tho fiu-rKy-Mt.iti; «odt*l, will bv 
Idt-ntlcal, but tlu* p.itlnt will bv dltfi-r«-nt. 

A mure reallHtl> dynumli. HHidvl tli.iii the two- 
Mt'ite one In ubiatned by iiu ludlnit the rota* 
tlon.ll dvnamic-a as a ■m.ill efie<t .ind allow- 
Init drag dependi nrv on lift. Iliia leads tc 
the Ktatv equations 


jY - v(e,h)lK(e.h) - ^,.'(e.h,L>) 

^1 
dt 

d7 . 1 


• v(e,hlsin 1 


(<tl) 


dl v(e,h) 


(i. - cos ■> ) 


where D|'(e,h,L) is the increment in drag- 

due-to-llft over thai. at L • mG. The control 
variable in (41) is L, and Che problem is 
nonslngular . This formulation of the problem 
is realistic enough for most purposes and was 
solved to first order in c by Ardema (1975, 
1976) by the method of matched asymptotic 
expansions. This solution was found to be In 
excellent agreement with a steepest descent 
solution, and it will be considered “exact" 
for the purposes of comparison. 

There have been two other analyses of the 
minimum t ime-to-cl imb problem by singular 
perturbation methods. Cnllse (1975) con- 
sidered the system 


^ - v(e,h)lF(e.h) - ^ Dj^’(e,h,D) 


flit) 

dh . 

dt 

V(e,h) 

sin y 

1 

(42) 

f,(c) 

& 


1 

(L - cos y) 



- 

dt 

v(e,h) 

i 
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Here, M Is lilt Hai li number whiih Is defined 
by M ■ v/a(h> wlicr4* i Is tbr spted of 
Hciind, I known fun' lion of altitude, f r 
formal. il Ions in whiili ) is not nwideied us a 
>late variable, t h< boundary tonJltlons on > 
do n<<t applv; a slmtlar 4t>mm«’nt applies to h. 
Ihe data describing the air* raft ire given in 
Ardtm.i (1975) and brvson, besni, and Hoff man 
(1969). Figure i shi>ws how the addlti'. e com- 
posite for altitude. h,(t). Is formed from 
its constituent fun. t ions a*.'ording to (IM). 
Ihis figure clearly .shows t h< stability of 
the boundary layur solutions and the nature 
of the composite solution: near t ■ 0, h.|(t) 

resembles hj(t): n*.ir t - t^ it resembles 
h|(tf - t): and awav from the boundaries it 

is identical to h,, It), t.* and are 

seen to be approximately i2 s and 20 s respec- 
tively. 

Table 1 compares the t lmes-ti>-cl imb as com- 
puted by the various methods. The additive 
composite solution (also energy-state) under- 
estimates the time by a substantial amount 
while the patched solution overest lm.it es it 
by an equally large anu’unt. A similar result 
was found for the more exact model (41) In 
Ardemn (1975). The good agreement between 
the patched and exact values lo of no signifi- 
cance since the patchcu solution Is an 
approximation of the two-state and not the 
exact solution. 


Table 1. Comparison of M i nlmum Tim e— to-Cl Imb 
Various Methods 


Tlme- 

to-cllmb 

8 


Energy-state (39) with (40) 105 

Two-state (29) with t ■ 1 130 


where 


11m fj(c) > 0 
e 


"Patched" approximation to two-state, 

(36) added to ends of solution of 

(34) with (35) 157 


1 Im 
e-»o 


f2(c) 

fl(t) 


0 


This formulation leads to multiple boundary 
layers. Breakwoll (1977) consideri'd the 
system (41) with e ■ 1 and the reclptrocal 
of maximum lift-to-drag ratio treated as a 
small parameter. This formulation Is appeal- 
ing hc-cause of the physical significance of 
Che small parameter. 


Tlie numerical example of Ardem.i (1975, 1976) 
will now be solved to Illustrate and compare 
the various solutions. 'Ihe boundary condf- 
C ions are 


M(0) ■ 0.5, h(0) ■ 12,200 meters, >(0) ■ 0 
M(tf) ■ 2.0, h(tf) - 24,400 meters, 7(tj) FREE 


Zeroth order additive composite 
approximation to two-state (38) 105 

Exact, first-order solution of (41) 
or steepest descent 162 


Figure 3 compares the paths in the (h,H) 
plane. All the paths agree qualitatively In 
that they have tlie same characteristic shape: 
1) an Initial steep dive of nearly constant 
energy to the region of the energy-climb path, 
during which velocity is gained at the expense 
of a loss in altitude; 2) an intermediate por- 
tion spent in the vicinity of the energy climb 
path (locus along which energy can be gained 
most quickly) where velocity Is accumulated 
with little change In altitude; and 3) a 
steep "soom-cllmb" of nearly constant energy, 
during which the excess speed Is traded for 


6 


.iltiiltdt gain. It in iliia cli ir.it t . r I at Ic 
'i t|<i t^liltli. In r.i' ' I mi»iH tliv pnihlt'n .iioi'n* 
•ihli- to ’• liitMi I .>r |>t'rtiir'' it Ion t«.‘c)in l<|ut‘H. 
iiii.inl It It Ivcly, ilio lwi»>t It. nolnllfn ip|K iri 
t.i give d i.preHint ii li'n .if ini- t-x.ii t 

..tiliit l>>n. T) tt /i rotii-.irdcr dtiilltlvi’ lonpo .iti. 
Ik n > better an .ipprox lm.it Ion to thv two-»t .ti 
Holiit Ion tli-in 1 h t)u* energv-Ht It I ( reUiict ti > 
solution. Again, .1 ilmll.ir rcHiilt has been 
fi'unJ lor tlie luiir* e'^n. t iniiJe I lArdfin.i, I**?!), 
riif ll•eflllne•a of tin /er.'t h-or.K r tolwi Ion 
In that it id rel.it ivelv eoNV to obt lin and 
I an be lilted .IN the banlH of a (lr»t-ordi r 
Nolutlm. Baaed on .\rdem.i'« (1975) reNultM. 
the flmt-ordt r Nolutlon nhould have very 
gi'od acc.irjt'V. 
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ALTITUDE. 1000 fwt 



Fig. 1. Sketch of Exact and Approximate Solutions for a Typical Case 


30.000 



Fig. 2. Additive Ciimpoaite for Altitude 
and Its Constituents 


ENERGY STATE. 
PATCHED APPROXIMATION 
TO TWO STATE 



Fig. 3. Comparison of Paths 
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16 AlHfracl 

A third order, nonlinear, singularly perturbed optimal control problem 1 


Is considered under assumptions which nsstire that the full problem is singu- 
lar and the reduced problem is nonsingular. The separation between the 
sinfjular arc of the full problem and the optimal control law of the reduced 
one, both of which are hypersurfaces in state space, is of the same order 
as the small parameter of the problem. Boundary layer solutions are con- 
str 'Cted which are stable and roach the outer solution In a finite time. 

A uniformly valid composite solution Is then formed fro i the reduced and 
boundary layer solutions. The value of the approximate soluti«in is that it 
Is relatively easy to obtain and does not Involve slnRular arcs. To illus- 
trate the utility of the results, the technique is used to obtain an 
approximate solution of a simplified version of the aircraft minimum t Ime- 
to-cllmb problem. numerical example is Included. 
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